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Abstract. Computer simulations have been employed in recent years to evaluate 
the configurational entropy changes in model glass- forming liquids. We consider two 
methods, both of which involve the calculation of the 'intra-basin' entropy as a means 
for obtaining the configurational entropy. The first method involves the evaluation of 
the intra-basin entropy from the vibrational frequencies of inherent structures, by 
making a harmonic approximation of the local potential energy topography. The 
second method employs simulations that confine the liquid within a localized region 
of configuration space by the imposition of constraints; apart from the choice of 
the constraints, no further assumptions are made. We compare the configurational 
entropies estimated for a model liquid (binary mixture of particles interacting via the 
Lennard-Jones potential) for a range of temperatures, at fixed density. 



PACS numbers: 
1. Introduction 

Whether a thermodynamic phase transition underlies the transformation of a 
supercooled liquid into an amorphous solid, or glass, at the laboratory glass transition 
temperature T g is among the central questions addressed by numerous researchers 
studying the supercooled liquid and glassy states. The notion of configurational 
entropy[|l], has played a significant role in attempts to define and understand the 
thermodynamic nature of the glass transition. In recent times, there have been various 
attempts to determine the configurational entropy of realistic liquids analytically and 
by computer simulations @, |, |, §, g, g [K| 0, 0, |E| 0, 0]J [L6], [T7j. The purpose 



of this paper is to compare two such methods that have been studied recently, namely 
the evaluation of the configurational entropy via the analysis of local potential energy 
minima or inherent structures (IS)|7|, g, ||, [H], [U], [16], |T7|], and by the calculation of 



basin free energies by confining the liquid within a localized region of configuration 
space by the imposition of constraints 0. These approaches, and results from their 
implementation, are described in the following sections. 

The model liquid studied is a binary mixture of 204 type A and 52 type B particles, 
interacting via the Lennard-Jones (LJ) potential, with parameters zab/^aa — 1-5, 
£bb/caa = 0.5, <Jab/o~aa = 0.8, and ubb/vaa = 0.88, and m^/m^ = 1, which has 
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been extensively studied as a model glass former [|18| , [T9| , |], [11], |T7||. Results presented 
in Sec. II are from molecular dynamics simulations, at a reduced density p = 1.2, 

lTf . Since the density is fixed, the 



which have been described in detail elsewhere M 



dependence on density is not always shown explicitly in the following. 



2. Configurational Entropy from Inherent Structures 

In the inherent structure approach 0], one considers the division of configurational space 
into basins of local potential energy minima. In practice such basins may be defined 
as the set of all points in configurational space that maps to the given local minimum 
under a specified local energy minimization procedure. Quite generally, one may then 
write the total partition function of the system as a sum of restricted partition function 
integrals over individual basins. Rewriting the partition function in this way introduces 
an entropy term associated with the number of local potential energy minima. With 
the expectation that configurations within a given basin are accessible to each other 
by thermal agitation while those belonging to distinct minima may not be, the number 
of distinct potential energy minima can be seen to be a measure of the number of 
physically distinct configurations or structures the system can adopt, i. e. a measure 
of the configurational entropy. 

Thus, the canonical partition function is re- written as a sum over all local potential 
energy minima, which introduces a distribution function for the number of minima at a 
given energy: 

Q N (p, T) = A"^— L-j / dr N exp (1) 
= ex P A ~ 3N I dr N exp (-/?($ - $ Q )) 

= J d<$> a exp(-(3($ a + Nf basin ($ a ,T))) 

= J d$ a exp (-/?($ Q + Nf basin ($*, T) - TS c (Q a ))) 

where <£> is the total potential energy of the system, a indexes individual inherent 
structures, $ Q is the potential energy at the minimum, V a is the basin of inherent 
structure a, fl(& a ) is the number density of inherent structures with energy $ a , and 
the configurational entropy S c = ks In Q (Note that here S c is a function of energy; the 
equilibrium average of this quantity displayed in Fig. 4 as a function of temperature). 

The probability of finding the system in the basin of an inherent structure of a 
given energy is given by the above as, 

P($ a , T) = 1 exp (-{l($ a + Nf basin (® a , T) - TS c {® a ))) . (2) 
Qn(P, T) 

The probability distribution P can be obtained from computer simulations, and 
offers a means of obtaining S c , provided one can estimate Qn (equivalently the free 
energy A(p, T) of the system) and the basin free energy f basing a, T). 
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The free energy at any desired temperature is obtained from thermodynamic 
integration of pressure and potential energy data from MD simulations ||. The 
absolute free energy A(p, T) of the system at density p at a reference temperature 
T r = 3.0 is first defined in terms of the ideal gas contribution A id (p,T) and the excess 
free energy A ex (p,T) obtained by integrating the pressure from simulations: 

A(p,T) = A id (p,T) + A ex (p,T), (3) 
PA id (p,T) = N(3 InA + lnp-l), 

cp dp (0 r P 
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Here, N is the number of particles, f3 = ksT, A is the de Broglie wavelength, and 
A° ex arises from the mixing entropy. A ex at a desired temperature may be evaluated by 
integrating the potential energy, E: 

(3A ex (p, (3) = (3A ex (p, p r ) + f E(p, (3')d(3' (4) 



Pr 

As observed in [|, |TTJ, the T dependence of E at the studied density is well described 
by the form E(p,T) ~ T 3 / 5 , in agreement with predictions for dense liquids [Til]]. A fit 
of the potential energy data to this form affords a means of extending with confidence 
the temperature dependence of E to T values where direct MD data is unavailable. 

The basin free energy fbasin(§a,T) is obtained by a restricted partition function 
sum over a given inherent structure basin, V a . For sufficiently low temperatures, one 
may expect the basin to be harmonic to a good approximation. In the harmonic 
approximation, we have 

3/3 1 3Ar " 3 

Pf=0 MtH + 7777 111 X i = P f therm + Pfvib, (5) 

where \ are eigenvalues of the Hessian or curvature matrix at the minimum. For 
individual minima, these eigen values are obtained by numerical diagonalization of 
the Hessian. The basin free energy can then be obtained either as a function of the 
inherent structure energy (by averaging free energies within individual energy bins) 
or as a function of temperature, by averaging all inherent structures sampled at a 
given temperature. (3f V ib is a slowly varying function of temperature (the temperature 
dependence is obtained by averaging over 1000, 100 inherent structures at T < 1., T > 1. 
respectively), and is fitted to the form j3f vi b(T) = fo + fi/T 2 which fits available data 
quite well. 

If the harmonic approximation to the basin free energy is accurate, inversion of Eq. 
(|), expressing S e ($ a ) in terms of P($ Q ,T), Q N (p,T) (or A(p,T)) and f ba sin{$a,T), 
for different temperatures T, should result in curves that overlap with each other, as 
S c (<& a ) is independent of T. Figure 1 shows the result of such inversion, which indicates 
that below T = 0.8, the various S c curves do overlap, while they do not at higher T. 
The procedure applied here is similar to, but improves upon, the procedure of shifting 




Figure 1. Configurational entropy as a function of inherent structure energy per 
particle, obtained from a range of temperatures. The overlap of curves for T < 0.8 
indicates that the harmonic approximation to the basin free energy is reasonable for 
T < 0.8. The solid line is a quadratic fit. 



unnormalized S c curves adopted in H, 10 . Thus, Fig. 1 indicates that a harmonic 



approximation to the basin free energy is not valid for temperatures higher than T = 0.8. 
The temperature dependence of the average inherent structure energy Eis, shown in 
Fig. 2, is consistent with this conclusion. As discussed in a simple expectation 
for the T-dependence of the average inherent structure energy in the harmonic regime 
is that Ejs ~ 1/T. Figure 2 shows that such a T-dependence is indeed valid at 
low temperatures, but breaks down for T > 0.8. However, this observation must be 
viewed in conjunction with two other observations about the topography of the inherent 
structure basins: (i) it has been demonstrated recently |2l| that the separation between 
'vibrational' and 'inter-basin' relaxation becomes reasonable for temperatures close to 
and below the mode coupling T c (~ 0.45 for the model liquid studied here), (ii) The 
difference in the potential energy of instantaneous configurations and the corresponding 
inherent structures is nearly linear with a slope of 3/2 for temperatures as high as 
T = 1.5, as shown in Fig. 3. Such a linear temperature dependence would normally be 
associated with harmonic behavior, which in the present case is misleading. 

The total entropy of the liquid S as well as the basin entropy Sb asin are evaluated 
as a function of density and temperature from the total and basin free energies. The 



Evaluation of configurational entropy of a model liquid from computer simulations 5 




_7 I _ l_ _^ _^ 

1 2 

T" 1 (Inverse Temperature) 



Figure 2. Average inherent structure energy vs. inverse temperature, showing that at 
low temperatures, the temperature dependence is well-described by T -1 while above 
T = 0.8 this dependence is not valid. 



configurational entropy S c (p, T) and the ideal glass transition Tia(p) are then given by, 

S c (p,T) = S(p,T) - S basin (p,T); S c (p,T IG (p)) = 0. (6) 

Figure 4 shows the configurational entropy so obtained as a function of T. 
By extrapolation, based on the assumption that the potential energy varied with 
temperature as T 3//5 , the ideal glass transition occurs at T = 0.294, in good agreement 



with estimates in IE, 11 



3. Constrained System Simulations 

An alternate approach to defining the basin entropy, which has been explored by 
Speedy||] is to impose constraints on a liquid to trap it in one of the basins it 
samples in equilibrium. A related approach has also been studied in [jnfl. With 
suitably chosen constraints, the calculated properties of the constrained system allow 
the evaluation of the basin entropy. A reasonable choice of constraint will restrict the 
system to a physically meaningful set of configurations related to each other without 
the need for configurational rearrangement. Further, such a constrained system should 
behave reversibly. In this work, the usefulness of one simple constraint is explored, by 
calculating the configurational entropy for a set of six temperatures at a fixed density of 
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Figure 3. Temperature dependence of the difference between the potential energy of 
liquid configurations and their corresponding inherent structures, which is linear with 
a slope of 3/2 upto very high temperatures. 



1.2, and compared with corresponding results from the inherent structure calculations 
described above. It is found that the constrained simulations result in comparable 
numbers for the configurational entropy from the inherent structure results. 

Ten sample configurations are chosen at each temperature, and the Voronoi 
tessellation is performed for each configuration. The Voronoi cell of each given particle, 
and the corresponding geometric neighbors, correspond to the cage a particle experiences 
at short and intermediate time scales. A configurational rearrangement of particles in 
the system will result in a restructuring of the Voronoi tessellation as well. Thus, the 
constraint of restricting particles to their Voronoi cells is an a priori reasonable choice. 
Hence, a constraint is imposed which confines each particle to its Voronoi cell during 
the Monte Carlo simulation from which the properties of this constrained system are 
evaluated. Each Monte Carlo simulation mentioned below is performed for 25, 000 Monte 
Carlo steps. The constrained system can be studied at any desired temperature; the 
temperature of the simulation from which the reference configurations are taken will 
be referred to as the Active temperature where there is need to distinguish these two 
temperatures. In order to estimate the configurational entropy, we must evaluate the free 
energy of the constrained system. This is done by thermodynamic integration |22, ^ 



from a reference system where each particle experiences a harmonic potential around 
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Figure 4. Configurational entropy vs. temperature, obtained from (a) Inherent 
structures (solid line), (b) constrained simulations, where the constraint is applied to 
equilibrated liquid configurations (filled circles, labeled 'Constrained Simulations 1'), 
and (c) constrained simulations, where the constraint is applied to inherent structures 
(filled squares, labelled 'Constrained Simulations 2'). Inset shows the extrapolation of 
S c which vanishes at T = 0.294. 



the initial configuration (Einstein cystal). Considering a potential energy function of 
the form, 

$(A, r N ) = (1 - X 2 )($ LJ (v N ) + $ c ) + A 2 C$> - r?) 2 (7) 

i 

where A is a tuning parameter that varies between and 1, is the Lennard- Jones 
potential of the unconstrained system, $ c is the constraining potential (which is zero if 
the constraint is obeyed and infinity if it is not), the corresponding free energy is given 
by 



A(X,p,T) = -k B T log 



A -3N J dr N exp(-f3$(\,r N )) 



(8) 



d\ 



The required free energy, A(X = 0, p, T) is related to that of the Einstein crystal (which 
may be calculated straightforwardly), by 

r 1 dA 
/o ~d\ 

where, from differentiating Eq. (§) with respect to A, 
dA 



A(X = 0,p,T)=A(X = l,p,T) 



-2A<$ LJ -C5>-r°) 2 >. 



(9) 



(10) 



Evaluation of configurational entropy of a model liquid from computer simulations 8 



The required average in the above equation is calculated by performing Monte Carlo 
simulations for a set of 11 A values. The values of dA/dX obtained are shown in Fig. 5. 
The free energy for A = 1 is 



A(X = l,p,T) = 3Nk B T 



logA — M^; 



All the above calculations are also performed using the inherent structures corresponding 
to the equilibrated liquid configurations mentioned above. With the free energy of 
the liquid evaluated as described in the previous section and the free energy of the 
constrained system obtained as described here, the configurational entropy of the system 
is given by 

Sc /k B = — (12) 

c/ B Nk B T Nk B T K ' 

where A cs is the free energy of the constrained system. The resulting configurational 

entropies are shown in Fig. 4. The S c values from the constrained system simulations 

are comparable with the inherent structure results, but the agreement is moderate. In 

particular, the constrained system results vary more weakly with temperature. 

To verify that the chosen constraint is a reasonable one, the free energies of the 

constrained system (for configurations from equilibrium runs at (Active) temperatures 

Tf = 0.629,0.484) are obtained independently at T = 0.1 and T = 1.2 from 

thermodynamic integration with respect to the Einstein crystal. Simulations are also 

performed for temperatures in between these two values, from which the temperature 

dependence of the potential energy is obtained. Using Eq. (f|), A(T = 0A)/Nk B T 

is calculated by integrating from T = 1.2. The difference of the directly calculated 

value and the one by integration is found to be —0.0197 for Tf = 0.484 and .0247 for 

Tf = 0.629. In other words, the constrained system appears to be reversible within 

the margin of error represented by these numbers. However, the discrepancy in the S c 

values between the constrained system and the inherent structure estimates is of the 

same order. Indeed, the discrepancy in the S c values for T = 0.484 and T = 0.629 is 

roughly the same as the discrepancy in the free energies above. It is likely that the 

sample of ten configurations used here is too small to obtain more accurate values. 



4. Conclusions 



Configurational entropy is obtained for a binary mixture liquid from analysis of 
inherent structures, and from estimation of the basin free energy via constrained 
system simulations. While the harmonic approximation used in the inherent structure 
approach to evaluate the basin free energy is in general questionable, the difficulty in 
the constrained system approach is the proper choice of constraint. The values for the 
configurational entropy obtained are comparable but show only moderate agreement. 
Further tests for the accuracy of the constrained system results, and more importantly, 
exploration of improved constraining methods are desirable for making a more stringent 
comparison of these two methods of calculating the configurational entropy. 
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Figure 5. Values of ^ for temperatures as marked. The filled symbols represent 
values for systems constrained with respect to the corresponding inherent structures. 
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